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Abstract 

Studying the substructure of jets has become a powerful tool for event discrimination and for 
studying QCD. Typically, jet substructure studies rely on Monte Carlo simulation for vetting their 
usefulness; however, when possible, it is also important to compute observables with analytic meth- 
ods. Here, we present a global next-to-leading-log resummation of the angular correlation function 
which measures the contribution to the mass of a jet from constituents that are within an angle R 
with respect to one another. For a scale-invariant jet, the angular correlation function should scale 
as a power of R. Deviations from this behavior can be traced to the breaking of scale invariance 
in QCD. To do the resummation, we use soft-collinear effective theory relying on the recent proof 
of factorization of jet observables at e'^e" colliders. Non-trivial requirements of factorization of 
the angular correlation function are discussed. The calculation is compared to Monte Carlo parton 
shower and next-to-leading order results. The different calculations are important in distinct phase 
space regions and exhibit that jets in QCD are, to very good approximation, scale invariant over 
a wide dynamical range. 
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I. INTRODUCTION 



The current success of the Large Hadron Colhder (LHC), its high center of mass energies, 
its significant dehvered integrated luminosity and its high-precision experiments has ushered 
in a new era of particle physics. Particles and jets with significant transverse boosts are now 
being copiously produced. An entire field of studying the substructure of highly boosted 
jets has grown up out of the study of these objects and many methods have been proposed 
to study QCD. In addition, procedures for discriminating QCD jets from jets initiated by 
heavy particle decays have been introduced and new measurements of these methods are 
being completed [HE]. To understand these methods in detail, most analyses have relied 
on Monte Carlo simulation as the basis of study. However, Monte Carlo simulations have 
limitations, and, where possible, it is vital to also compute the observables to higher orders 
in QCD so as to have another handle on their behavior. 

An important contribution to this effort of computing jet observables is resummation 
of large logarithms that arise in fixed-order perturbation theory. Jets are objects that 
are typically dominated by soft or coUinear emissions and so it is necessary to resum the 
logarithms that exist for an accurate prediction of an observable. Very recently, groups 
have computed resummed contributions to light jet masses at hadron colliders [3] and N- 
subjettiness [H [5] in color-singlet jets at the LHC [6]. Ref. [6] in particular relied on the 
factorization of color singlet processes at hadron colliders to reinterpret results from e^e~ 
colliders. Computing the resummed contribution to generic observables at hadron colliders is 
made more difficult by the color flow throughout the collision which can destroy factorization. 
To avoid discussion of these issues, here we will only consider jet observables at e^e~ colliders. 
In this paper, we will discuss the resummation of the angular correlation function introduced 
in [7] using soft coUinear effective theory (SCET) [8tiTT]. 

The angular correlation function Q{R) was defined in [7] as 



for studying the substructure of jets at the LHC. Ai?jj is the boost-invariant angle between 
particles i and j, the sum runs over all constituents of a jet and is the Heaviside theta 
function. The angular correlation function has distinct properties for scaleless jets versus jets 
with at least one heavy mass scale. In particular, any structure in the angular correlation 
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function should be distributed roughly as i?^, where D is a constant, for a scaleless jet. 
It was shown that by exploiting the different behavior of the angular distribution of hard 
structure in QCD jets versus jets initiated by heavy particle decay, an efficient tagging 
algorithm could be defined. 

Ref. [12] continued studying the properties of the angular correlation function, focusing 
on average properties of QCD jets. It was shown through simple calculations that, for QCD 
jets, the angular correlation function averaged over an ensemble of jets should approximately 
scale as 

{g{R)) ^ , (2) 
where the angle brackets are defined by 

Deviations from are due to the running coupling and higher order effects. The introduc- 
tion of an ensemble averaged angular correlation function allows for a rigorous definition of 
the dimension of a QCD jet which is also infrared and coUinear (IRC) safe. This dimension 
is defined to be the average angular structure function (A^) and is the power to which the 
average angular correlation function scales with R: 

For QCD jets, (A^) ~ 2. In [T2j, it was also shown that the scaling of non-perturbative 
physics in R is distinctly different, and this was used to determine the average energy density 
of the underlying event. 

Here, we will continue the work of [12] and compute the average angular structure function 
by resummation within the context of SCET. Our analysis is only truly appropriate at e^e~ 
colliders, but we expect that the largest effect in going to hadron colliders is the contribution 
of the underlying event. For this calculation, we introduce generalized correlation functions 
Qa{R) parametrized by an index a: 

The form of the angular correlation function is similar to jet angularity [131 E] ^-^id close 
in form and spirit to an event shape introduced in [? ]. However, for our purposes here, we 



choose to index the parameter a such that the angular correlation function is IRC safe for 
all a > 0. In the small angle limit, this reduces to Eq. [ijwith a = 2 (up to normalization). 
The parameter a allows for a study of the behavior of the angular correlation function with 
angular scales weighted differently. Analogously to the angular structure function, we define 
a generalized average angular structure function 

^^^-^ = -dl^ ■ 

The calculation and interpretation of average angular structure function will be the focus of 
this paper. 

In Sec. [IT], we discuss the factorization of jet observables in SCET and the computation 
of the angular correlation function including global next-to-leading-log (NLL) contributions 
for jets defined by a kj-tjpe algorithm [TB] . The existence of a factorization theorem for 
the angular correlation function is non-trivial. We will discuss the consistency conditions 
that the angular correlation function satisfies for factorization. We will also briefly discuss 
how the results obtained here can be used in a calculation of the angular correlation function 



at the LHC In Sec. Ill we compare the SCET calculation to a next-to-leading-order (NLO) 
calculation of the angular correlation function. Resummation and fixed-order corrections 
affect different parts of distributions and so the differences between the resummed calculation 
and the fixed-order result give some sense as to the importance of these effects. This analysis 



leads to Sec. [TVj were we present a comparison between the SCET calculation and the output 
of parton shower Monte Carlo. We observe significant differences between SCET and Monte 
Carlo, but higher fixed order effects are substantial. We discuss some of the uncertainties in 
the parton shower studying the effect of the evolution variable on the value of the angular 
structure function. Finally, we present our conclusions in Sec. |V] 

II. SCET CALCULATION 

SCET is an effective theory of QCD in which all modes of QCD are integrated out except 
those corresponding to soft or coUinear modes. CoUinear and soft modes are defined by their 
scaling with power counting parameter A: 

coUinear ~ (A^, 1, A) , 
soft ~ (A^A^A2) , 



which is the scahng of the +, — and transverse components of the momenta, respectively. 
A is a parameter that is defined for a particular process or observable; for example, for 
computing the distribution of jet masses, A ~ rajj-p^j ^ 1. The fact that A ^ 1 allows for 
a systematic expansion in powers of A. Higher order terms in A are power suppressed (much 
like the subleading terms in the twist expansion). 

For an event shape observable O that factorizes, the cross section can be written in the 
schematic form: 

da 



s{o-^), (7) 



where H{^) is the hard function, which matches the full QCD result at a scale /x, J(/i; n^, O) 
is the jet function for the contribution to the observable O from nj-collinear modes and 
S{^\ O) is the soft function for the contribution to the observable O from the soft modes. 
® represents a convolution between the jet and soft functions. All functions depend on the 
factorization scale /i. 

Factorization of jet observables in SCET was first exhibited in [TTIIIH]. Ref. [18] computed 
individual jet angularities to NLL in e^e~ collisions. It was shown that factorization of the 
cross section for jet observables in e~^e~ — )■ N jets has the form 



= . . . ,nAr;/i) 



dOi--- dO 
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N 

Sn,...nAOl,...,OM;fi) n -^(Z^) '(8) 

j=M+l 



where M < of the jet observables Oi have been measured. Jet directions are denoted by 
rii and J{Oi] fi) is the jet function for a jet in which the observable Oi has been measured 
and J(/i) is a jet function for a jet which has not been measured. We will refer to these as 
the measured and unmeasured jet functions, respectively. A similar nomenclature will be 
used for the soft functions. Jet algorithm dependence and jet energies have been suppressed. 
An important point from [TH] is that factorization requires that the jets be well-separated; 
namely, that 

tan^ 

= wF ' • 

where ipij the is angle between any pair of jets i, j and Rq is the jet algorithm radius. We will 
assume that this condition is met in the following and leave any discussion of subtleties to 
|18j . A non-trivial requirement of the factorization is the independence of the cross section 
on the factorization scale /x. This requirement leads to a constraint that the sum of the 



anomalous dimensions of the hard, jet and soft functions is zero. We will show that this 
holds for the angular correlation function. 

We will use the results of the factorization theorem proven in [T71 |TS| to compute the 
distribution of the angular correlation function from Eq. |5| In particular, we are interested 
in the ensemble average of the angular structure function as defined in Eq. [6] Note that this 
observable is independent of any normalization factor of the angular correlation function; 
thus, with the goal of computing the average angular structure function, it is consistent to 
ignore factors that are independent of Qa and the angular resolution parameter R. Thus, 
for the purposes of this paper, we can ignore the overall factors in the factorized form of 
the cross section of the hard function and the unmeasured jet functions. In this case, the 
factorized form of the cross section becomes 

{Gal, ■ • ■ , QaM] /i) , (10) 

where C(/x) is independent of Qa and the resolution parameter R. 

In this section, we present a calculation of the jet and soft functions for the angular 
correlation function for jets defined by a algorithm. We first argue that the angular 
correlation function is computable in SCET and relate its form at NLO to the form of jet 
angularity at NLO. This comparison will allow us to relate the calculation of the angular 
correlation function to the work in [18J. We then present a calculation of the measured jet 
and soft functions of the angular correlation function. From these results, we can determine 
the anomalous dimensions of the jet and soft functions and will show the consistency of the 
factorization relies on a non-trivial cancellation of dependence on the angular resolution R 
between the jet and soft functions. We can then resum up to the next-to-leading logs of 
the jet and soft functions by the renormalization group. Note that we do not attempt to 
resum non-global logs [19] that arise due to the non-trivial phase space constraints of the 
jet algorithm or the angular correlation function. From the resummed expression of the 
angular correlation function, we find the ensemble average and compute the average angular 
structure function numerically. 

It should be stressed that non-global logarithms are ignored in this study. The angular 
correlation function for a jet requires several phase space constraints; the jet algorithm, soft 
jet vetoes, the resolution parameter R, etc. These provide numerous sources for non-global 
logarithms which cannot be resummed analytically. The study of non-global logarithms in 
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QCD cross sections is a subtle and evolving story. For recent work in this direction, especially 
in the context of non-global logarithms from jet clustering see, for example, [2nH2S]- It is 
outside the scope of this paper to discuss non-global logarithms further. 

A. Factorization of the Angular Correlation Function 

Factorization of jet observables requires that soft modes only resolve the entire jet and 
not individual coUinear modes contributing to the jet. Angularity Ta is a one-parameter 
family of observables defined as ^131 [H] 

where J is the jet, p±i is the momentum of particle i transverse to the jet axis and rji is the 
rapidity of particle i with respect to the jet axis: 

6- 

rji = - log tan . (12) 

Angularity is IRC safe for a < 2. The separation of soft and collinear modes in angularity 
is simple to show. To leading power in A, 

= + r! , (13) 

where C and S represent the collinear and soft modes, respectively. Note that the soft 
modes do not affect the location of the jet center to leading power in A. Factorization of 
angularities exists only for a < 1 due to the presence of logarithms of rapidity; however, 
recently it was shown that these logarithms can be controlled |26l |2Z|- We will show that 
angularity and the angular correlation function have similarities which will allow us to use 
many of the results from [18] here. 

To justify the use of SCET for computing the angular correlation function, we must first 
show that the angular correlation function does not mix soft and collinear modes. This 
argument was presented in [12j (based on arguments from [28]), but we present it here for 
completeness. In terms of soft and collinear modes, the angular correlation function can be 
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expressed as 

GaiR) = ^ E ^^^^ % Y©(^ - 



Note that, to NLO, there is no soft-soft correlation contribution to the angular correlation 
function because such a term would require the radiation of two soft gluons which first occurs 
at NNLO. To accuracy of the leading power in A, we can exchange the coUinear modes with 
the jet itself in the coUinear-soft term. Explicitly, 

Ocs = Ojs + 0{X) , (15) 

as the angle of the soft modes with respect to the jet center scales as 9js ~ 1- Appropriate 
for NLO or NLL, the angular correlation function can be written as 



G^iR) = ^ E ^^^^ Y'^iR- %) 

+ ^J2^ssmejstan--'^-fe{R-9js) . (16) 

s 

Thus, the collinear and soft modes are decoupled to leading power and so the angular 
correlation function is factorizable, and hence computable, in SCET. 

To NLO, a jet is composed of at most two particles, so the form of many observables 
simplifies substantially at this order. The form of the angular correlation function from Eq. [5] 
was chosen so as to be similar in form to angularity. The contribution to the angularity and 
the angular correlation function from collinear modes is distinct. The measured jet functions 
will need to be recomputed for the angular correlation function. However, the contributions 
to the angularity and the angular correlation function from soft modes are simply related: 

^^f (^) = 7^ sin^jstan^-i ^e(i? - 9js) = rf_„e(i? - Ojs) • (17) 

This observation will allow us to recycle the soft function calculation for angularity for the 
angular correlation function. 
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An important point to note here is that the scahng of the angle between coUinear modes 
i and j goes hke 9ij ~ A. Thus, to leading power, the angular correlation function for the 
coUinear-coUinear contribution can be written as 

= ^ E ^^^^ % y ® - %) 

= ^2 E ^^^^ t^"'' y - ^^^■) ■ (18) 

We will use this form of the coUinear-coUinear contribution to the angular correlation func- 
tion for computing the measured jet functions. 

Throughout this paper, we will only consider jets with a single coUinear sector. Small 
values of the angular correlation function are not enough to guarantee that the jet has 
only a single coUinear sector, however, we believe that contributions from multiple coUinear 
sectors is subdominant. Our reasoning is as follows. First, at large values of i? ~ i?o, the 
angular correlation function is essentially an angular- weighted jet mass measure. In this 
case, additional coUinear sectors would be correlated increasing the value of the angular 
correlation function substantially. At small R, there are two options: either coUinear sectors 
are still correlated or they are uncorrelated. If the coUinear sectors are still correlated at 
small R (they are within an angle R of one another), then logarithms of this angle will 
appear. However, the logarithms of the angle between the separate coUinear modes should 
be subdominant to the logarithms of the resolution parameter of the angular correlation 
function, R. We do not attempt to rcsum the latter logarithms in this paper. If the coUinear 
sectors are uncorrelated and separated by an angle larger than R, then the jet effectively 
breaks up into several smaller jets each with similar scaling properties as i? ^ 0. In sum, 
we expect the effect of additional coUinear sectors to be significantly subdominant so as to 
be consistently ignored in this study. We beheve the absence of fixed-order terms is a more 
important omission. 

B. Measured Jet Functions 

The leading power contribution to the measured jet functions at NLO comes from two 
coUinear particles which are clustered in the jet and can be computed from cutting one- 
loop SCET diagrams. The phase space integrals can be extended over the entire range of 

9 



momentum for the coUinear particles in the jet as long as the contribution from the zero 
momentum bin is subtracted [22] • In particular, we consider a jet with light cone momentum 
/ = 0) which splits to two collinear particles with light cone momenta q = , q_|_) 

and I — q = {l^ — q^, u) — q^ , — q±)- The zero-bin subtraction term can be determined from 
the measured jet function by taking the scaling q ~ A^. We will refer to contribution to the 
jet function that does not include the zero-bin subtraction as the naiVe contribution. 

To compute the measured jet function, we will need to enforce phase space cuts from 
the jet algorithm and the observable. We will compute the jet function for a kx-type jet 
algorithm as defined by a jet radius Rq. At NLO, all algorithms are the same and two 
particles are clustered in the jet if their angular separation is less than Rq. This leads to 
the phase space constraint 

e., = 9 I cosi^o - , q • " = 9 ( tan^ ^ - J^""' ^ , (19) 



|q|v/(l-q)V V 2 q-(uJ-q-)\ 

where on the right, the leading scaling behavior was kept. The jet algorithm constraint for 
the zero-bin subtraction term is then 

eS = e (tan= f - ^) . (20) 

The phase space constraints for the angular correlation function are more subtle. The 
^-function which constrains a jet to have angular correlation function Qa, = S{Qa — Qa), 
is 

Sr = s(ga- ^-\u. - q-f-{q-Y-'\q^r"Q (tan^ f - -^^^j) ^ (21) 

where R is the resolution parameter of the angular correlation function. For a fcy-type jet at 
NLO, the angular correlation function vanishes if i? > Rq\ thus, we will assume that R < Rq 
in the following. This 5-function can be decomposed depending on the value of 9-function 

as 

Sr = 6(g^- u^-'ico - q'r-{q~Y~-'\q^r''Q ftan^ f - f""' \ 
V V 2 q-{(^-q-yj 

= 5{g^- u-\uj - q-f-^q-f-^'^q-'r'') 9 ftan^ ^ 



q {u-q Y 



+ ^(g.)ef f""' ■ (22) 

\q [uj-q f 2 J 

The 5-function for the zero-bin subtraction term is found by taking q 
6f =6{G^- u-\q-f-l\q^r/') 9 (tan^ f " ^) + (^") ® " f ) ' (23) 
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FIG. 1. SCET Feynman diagrams contributing to the quark jet function. 



1. Measured Quark Jet Function 

The naive contribution to the measured quark jet function can be computed in dimen- 
sional regularization from the diagrams shown in Fig. [TJ 



X 



2n J (27rY V g ^ tu - q- 



2ti5 (q+q- - ql) Q(q-)Q(q+)2n5 ( - q+ ^] 

V ^-Q ) 



xe(a; - g-)e(/+ - g+)e f tan^ ^ 



q^u? 



q [u-q- 



X 



6 {g^ - u;--\u; - q-f-^iq-f-^'^q-'T") Q tan^ ^ - 



R q+u'^ 



2 q [00 — q 



-\2 



q+UJ^ 2 R 



q [u-q 



(24) 



We take d = 4 — 2e. The coefficient to the 6{Qa) term can be found by integrating over Qa- 
The terms that remain are +-distributions, which integrate to zero. The zero-bin subtraction 
term follows from taking the scaling limit g ~ of the naive jet function above: 



dl+ 1 f d'^q 1 



2ti /^ 



{2'nYq 

x2n5 (/+ - g+) 0(/+ - g+)0 ( tan^ 



-27rS {q+q- -ql)Q{q-)e{q' 



6 {g^ - co-\q-y-^/\q+Y^^) e tan^ 



2 q- 

R g^ 



tan 



R 



(25) 



The term proportional to S{Qa) is scaleless and integrates to zero in pure dimensional regu- 
lation. 

Employing a MS scheme, we find the measured quark jet function for /c^-type jet algo- 
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FIG. 2. SCET Feynman diagrams contributing to the gluon jet function. Diagrams (F) and (G) 
have mirrored counterparts which are not shown. 



rithms of Qa to be 



2lT 



a I 31 a log 

\ \ ^ 

- 1 e2 2 e a - 1 e 



a - 

1 /e(a 



1 tan^ f 



a 



le 




where J^^iGa,^*^) consists of terms that are finite as e — )■ 0. These terms are presented in 
Appendix |Aj The definition of the +-distribution is also given in Appendix |Xj Note that 
the 1/e terms for the angular correlation function are the same as those for angularity from 
with a — 7- 2 — a plus an additional term of the logarithm of the ratio of scales; the 
resolution scale R and the jet radius Rq. This term contributes to the anomalous dimension 
of the jet function. In principle, these logarithms could be attempted to be resummed. 
However, note that the resolution scale R can never practically be parametrically smaller 
than the jet radius Rq, so these logarithms never become large. Thus, we will not worry 
about resumming these logarithms. 

12 



2. Measured Gluon Jet Function 



The naive contribution to the measured gluon jet function can be computed from the 
diagrams shown in Fig. [2] 

2 q+q 



X <^ UfTr 1 



1 — e ut 



q^q- 



x0(g-)e(g+)e(u; - g")e(/+ - g+)0 (^tan ^ 



oo — q^ 
2R0 q'^ 



ujI+ 
2 



X 



6 {Go. - UJ--\co - q-f--{q-f-'^l\q-^r'') tan^ 



q [uj-q 
R 



-\2 



q^u 



q~{uj - q-y 



q+io^ 2 R 
^-tan 2 



(27) 



The coefficient of the 6{Qa) term can be found by integrating over Qa- The terms that 
remain are +-distributions, which integrate to zero. The zero-bin subtraction term follows 
from taking the scaling limit / — g ~ g ~ of the naive jet function above: 

dl+ 1 f d'^q 1 



271 /^ 



(27r)'^g 



-27r5 (gV-9l)e(g-)e(g^ 



X 



2n5 (/+ - g+) 0(/+ - g+)0 tan 



2 ^0 



X 



5(^.-c.-^(g-)^-'^/^(g+r/^)e t 



2 g- 

,R g^ 



;an 



q~ 



+ (5(6;«)0( ^-tan^l 



(28) 



The term proportional to S{Qa) integrates to zero in pure dimensional regulation. This 
zero-bin subtraction term is exactly the same up to color factors as the quark jet function 
zero-bin subtraction. 

Employing a MS scheme, we find the measured gluon jet function for the kx-type jet 
algorithms of Qa to be 



271 

+ 



a — le^ 2t a — I t 



Ca - tan 



log 



tan^ 



e{a - 1) V 



(29) 
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where J^{Qai^^) consists of terms that are finite as e — )■ 0. These terms are presented in 
Appendix |X} Pq is the coefficient of the one-loop /3-function: 

/3o = yC^A - \nf , (30) 

with Tij = |. As with the quark jet function, the 1/e terms are the same as those for 
angularity from [18j with a ^ 2 — a plus an additional term of the logarithm of the ratio of 
the resolution parameter R to the jet radius Rq. 



C. Measured Soft Function 



As shown above, there is a simple relationship between the form of angularity for soft 
modes and the angular correlation for soft modes. This relationship will allow us to use 
the results from [TH] in computing the measured soft function for the angular correlation 
function. First, we consider the phase space constraints from the jet algorithm and the 
angular correlation function. For the kr jet algorithm, soft radiation must be within the jet 
radius Rq of the jet axis to be included: 

= (tan^ y - ^) • (31) 

The ^-function that constrains the soft modes to contribute an amount Qa to the angular 
correlation function is 



5 



5 (g^ - cu-\k-Y-''/\k+r/^e (^tan^ f - ^) ) 
= 6{g^- u-\k-f'-/\k^r'') e (tan^ f - ^) + (^-) - tan^ f ) -(32) 
The measured soft function of a gluon emitted from lines i and j into a jet is 



s [g^ - uj-\k-y-^/\k+r/') e (^tan^ f - ^) 



+ 5(^«)e(^-tan2^ 



(33) 
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Note that the integral proportional to S{Qa) is scaleless and so vanishes in pure dimensional 
regularization. Also, the integral is only non-zero if R < Rq and so the 0-function from the 
jet algorithm is redundant. Thus, we can write the soft function as 



omeas 



{Go) 



(34) 



This is the same form of the measured soft function as for angularity with a jet radius equal 
to R which was computed in [TH] . Up to terms that are suppressed by from Eq. [oj the 
measured soft function for jet i is 



Cmeas ( f 

o [yi, 



271 ' a - 1 
-2 



1 1, /i^ tan2(°-i) I TT 



00^ 



2 1 „/i2tan2(--i) 
12 + 2^°^ 



1 , , /iHan2(-i)f\ e(^. 

- + log 7^2—2 



SiGia) 

(35) 



where is the square of the color in the jet. 



D. Anomalous Dimensions and Consistency Conditions 

A non-trivial requirement of the factorization is that the physical cross section should 
be independent of the factorization scale /i. A consequence of this is that the anomalous 
dimensions of the hard, jet and soft functions must sum to 0. The requirement is 

0= (7H(/i)+7r""^^^(/^)+ E iM)siGa)+ ^jM-, f^) + ir^Gi,-, f^)) , (36) 

\ i^moas / iGmeas 

where jh, 7s and 7j are the anomalous dimensions of the hard, soft and jet functions. The 
fi dependence must be summed over the measured and unmeasured jet and soft functions. 
The sum of the hard, unmeasured soft and unmeasured jet anomalous dimensions to NLO 
is 

i^meas lemeas * 2, iSmeas 

where 7, depends on the flavor of the jet: 

^^=2^^-' ^^ = V 6 = 2^^°' (^^^ 
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for quark and gluon jets respectively. We will show that the measured jet and soft function 
anomalous dimensions for the angular correlation function are exactly what is required to 



satisfy Eq. 36 



The anomalous dimensions of the measured jet or soft functions are given by the coefficient 



of the 1/e terms from Eqs. 26, 29 and 35 The anomalous dimensions of the quark and gluon 



jet functions can be written collectively as 



a.. 



a 



fi^ tan^ f 

log^ + log- — ^ 
a - 1 uf tan^ ^ 



+ 7* 



TT a — 



(39) 



where 7, is defined in Eq. |38j Note the non-trivial dependence of the anomalous dimension 
on both the jet radius and the resolution parameter of the angular correlation function. The 
anomalous dimension of the measured soft function for a quark or gluon jet is 



75 (t'aJ 



^rp2 

TT a 



;x2tan2("-i)f 

6{g^)iog^ , — ^-2 



(40) 



As mentioned earlier, jet angularity is not factorizable for a = 1 and here we see that the 
anomalous dimensions of the angular correlation jet and soft functions become meaningless 
for a = 1, signaling a breakdown of factorization. For the angular correlation function, we 
are most interested in a = 2, so we will not consider this issue further here. 

Summing over the measured jet and soft function anomalous dimensions, we find 

,2 



E (7J.(ei;A.)+7r"(ei;*')) 



iSmeas 



iSmeas 



uf tan^ ^ . 



}^ 7. s{g^) (41) 



Note that there is a non-trivial cancellation of the angular correlation function resolution 
parameter R between the jet and soft functions. This contribution exactly cancels that from 
the hard and unmeasured jet and soft functions in Eq. |37| consistent with the factorization 
requirement. 



E. Resummation and Averaging 

To proceed with the resummation to NLL of the jet and soft functions, we will make a 
few observations. First, as mentioned earlier, because we are ultimately interested in the 
average angular structure function, we can ignore factors in the resummed cross section that 
are independent of Qa or the resolution parameter R. Thus, we will not discuss nor resum the 
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hard function nor the unmeasured jet and soft functions. Also, we will only consider a single 
measured jet in an event. This prevents a study of inter-jet correlations of the angular 
correlation function, but for this paper we are most interested in the intra-jet dynamics. 
Anyway, the existence of factorization of jet observables essentially trivializes correlations 
between jets since it implies that correlations can only come from the soft function. From 
these observations, we only need to resum the measured jet and soft functions of a single 
jet. 

With these considerations, we will need to compute the convolution between the measured 
jet and soft functions: 

^ oc I J{g^ - Q'^- /X,, ix)S{Q'^- ixs. /^) , (42) 

where /ij and [is are the jet and soft scales respectively. We refer the reader to [IB] for 
the details of generic NLL-level resummation. Here, we will use the results collected there 
appropriate for the angular correlation function. The resummed differential cross section for 
the angular correlation function of a single measured jet at NLL is 



dCfa \ U / \ U I 1 ( — — Us) 

(43) 

u is the — component of the jet's momentum and je is the Euler-Mascheroni constant. The 
functions Uj, Us, Kj, Kg, fj and fs are written in detail in Appendix [Bj They depend on 
the jet and soft scales and the factorization scale fi. The jet and soft scales will be in general 
sensitive to the value of Qa and the resolution parameter R. At very small values of Qa, the 
resummed distribution can become negative and in general will need to be matched onto 
a non-perturbative shape function in that region. We do not attempt to correct the shape 
at very small Qa and instead just set the cross section to zero where it would otherwise be 
negative. 

The average angular correlation function can then be computed from the cross section in 



Eq. 43 by integrating over Q^- 
/•^"" da 

{Ga{R)) = / dg^ —g, 

Jo 



T{-uj - us) g'^''-''-' 
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tan*^ — 

where Q^^^ = — is the maximum value of the angular correlation function for a jet with 
two constituents. We choose the scales /xj and /is so as to eliminate the logarithms that 
remain in the resummed distribution. The choice of these scales can be seen from the form 
of the fj and fs terms as given in the appendix. We find 

= -QT , = ^^^^ . (45) 

With this choice of scales, the average angular correlation function simplifies: 

{GaiR)) oc / dg^ — ^ [l + fj + fs] . (46) 

Jo J- 1— i^j — ^s) 

Note, however, that there is non-trivial dependence on Qa in the functions uj, us, Kj and 

Ks- Finally, to determine (AQa), we compute 

<-^")-^- (-) 

Plots of the average angular structure function as computed in SCET and compared to 
Monte Carlo and NLO corrections will be presented in the following sections. 

While resummation of the angular correlation function is necessary for an accurate de- 
scription of the singular regions of phase space, it is not obvious how important resummation 
is for the average angular correlation function. Resummation of the distribution da/dQa 
tames the singularity at small values of Qa and produces a peak. The resummation contri- 
bution to the distribution is most important near the peak while fixed-order contributions 
are most important in the tail, at large values of Qa- However, the average angular corre- 
lation function is sensitive to both the resummed and fixed-order contributions. To get a 
sense of the importance of the resummed contribution, we can compare the location of the 
peak in da/dQa to the maximum Qa value possible for a jet with two constituents. If the 
ratio of the location of the peak to the maximum value of Qa is small, resummation effects 
are minimal while if that ratio is large, then the resummed contribution dominates. 

The comparison of the location of the peak in J?- to the maximum value of Qa is shown 
in Fig. [3j Here, we have set a = 2 for illustration and the plot shows how the location of 
the peak relative to the maximum value changes as the resolution parameter R decreases. 
When R = Ro, the angular correlation function is just the jet mass and the ratio is relatively 
small for both quark and gluon jets. However, as R is decreased from Ro, the ratio increases, 
reflecting the greater importance of the resummed contribution with respect to fixed-order 
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FIG. 3. Plots of the ratio between the location of the peak in to the majximiini VctluG of Qa 
over a range in R. For illustration, a = 2 and the red (blue) curve is quark (gluon) jets. The jet 
radius is i?o = 1-0 and we have set the hard, jet and soft scales as in Eq. |45} 

corrections. Thus, we expect that the fixed-order contribution to the average angular corre- 
lation and structure functions is largest at i? ~ -Rq while the resummed contribution becomes 
more important at smaller R. This will be discussed later when comparing the resummed 
calculation to Monte Carlo and NLO calculation of the average angular structure function. 




1. Low est- Order Expansion 



Before continuing, it is illuminating to expand the angular correlation function to lowest 



order in the coupling as- To do this, we will need to expand Eq. 43 to 0{as)- The form of all 
of the functions in Eq. 43 are given in Appendix |B] and, in particular, the expansions of the 
Gamma, harmonic number and polygamma functions are needed. The necessary expansions 
are given in the appendix. To leading order in a^, we find 



da 



oc 



271 ' 



2 R 



^ 4 , ^ 11 , tan „ 
4 log tan — log ^Q, I Cj + lo,- 



a 



a 



tan 



2 -Ro 



(48) 



where the factor q depends on the flavor of the jet: 



3 

2 ' 



/3o 



2Ca 



(49) 



To compute this, we have set the jet and soft scales so as to minimize the logarithms that 



appear in the cross section as defined in Eq. 45 In this expression, note that the non-cusp 
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piece of the anomalous dimension of the measured jet functions appears in the term in 
parentheses. 

From this expression for the cross section differential in the angular correlation function, 
we integrate over Qa to compute the average angular correlation function. To 0{as) we find 



„a R 



OC 



a R 
t X 







4 log tan 



R 



a 



log Ga 



tan^ f 



27r ' 



1 + log 4 



a 



Ci_ 

4 



1 tan^ I 



+ 0{a 



(50) 



Any overall factor independent of R does not affect the average angular structure function 
because 



{AGaiR)) 



d\0g{Ga) R d{Ga) 



(51) 



dlogR {Ga) dR 

To lowest order, the average angular structure function is independent of as and only de 



pendent on the color of the jet through the q term. Eq. [50] results in the average angular 
structure function of 



(AGaiR)) 



R 



sin R 



a 



4 + 4 log 4 



Ci + log —-^ 



(52) 



Eq. [52] contains much of the physics that we expect affects the form of the angular 

contains an 



structure function. The naive expectation for {AG a) is {AG a) — ct. Eq. 52 



(9(1) correction to this result that is negative. This was interpreted in [12] as an effect 
due to the running coupling. However, here, this is probably not the source of this effect 
because even for fixed coupling the negative term exists. This is instead probably due to 
SCET itself because only collinear and soft emissions are included with respect to full QCD. 
Including all terms in the resummed result should decrease the average angular structure 
function further due to both the running coupling and because an arbitrary number of soft 
and collinear emissions are considered. 

Also, note that the term Cj is larger for quarks than for gluons with sufficiently many 
flavors of quarks: 



3 ^ 11 Nf 
2 - "e" ~ "9~ 



--9 ' 



(53) 
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for Nf > 3. This implies that, for sufficiently many flavors, {AQa)g > {AGa)q, an observation 
that was also made in PI2]. There, this was attributed to the fact that gluons have more 
color than quarks and so radiate more at larger angles, effectively decreasing the strength of 
the collinear singularity with respect to quarks. We expect that the resummation magnifles 
the distinction between quarks and gluons. 

Another interesting observation to be made about the form of the angular structure 
function is that, to this order, it is Lorentz invariant. We then expect that all jets, regardless 
of energy (so long as it is above the hadronization scale of QCD), have an angular structure 



function that deviates only slightly from the form in Eq. 52 In particular, note that Eq. 52 
is the inflnite jet energy limit of the (all-orders) angular structure function. The contribution 
of higher orders to the angular structure function would contain prefactors of as{^) which 
would vanish as — ?■ oo. If we ignore the flnite R terms from the expansion of sine and 
tangent, {AQa{R)) is very flat, signifying very near scale invariance over a large dynamical 
range R. Flatness is only broken by a term that goes like 1 / log which is only important 
at very small R/Rq. 

It is accurate, then, to represent the angular structure function in the form (again, ig- 
noring the flnite R terms from sine and tangent) 



(Ae?„(i?)) ^«-7ASF , (54) 



where 7asf might be called the anomalous dimension of a QCD jet and is independent of 
a. This anomalous dimension is a robust quantity that is intrinsic to the flavor of the jet 
and properties of QCD. Measuring this property of the angular structure function in data 
would be very interesting. It is important to note, however, for all of the above comments, 
0{as) contributions to the average angular structure function have been ignored. These are 



expected to be comparable in size to the second term in Eq. 52 Note in particular that 
NNLO contributions can be just as, or even more, important than the contributions from 
resummation. Indeed, for jets with three constituents, it was computed in [12] that the 
effect at this order is to increase the average angular structure function. 
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F. Non-Perturbative Physics Effects 



In addition to the perturbative physics contribution to the angular structure function, we 
would also like to understand the effects from non-perturbative physics. For jets produced 
in an e^e~ collider, the dominant non-perturbative effect is from hadronization. A simple 
physical argument can be used to determine how hadronization affects the angular structure 
function. The partons created from the parton shower will be connected to one another 
by color strings which stretch across the event. After the termination of the parton shower 
at an energy scale of about 1 GeV, these color strings are allowed to break to create a 
quark-antiquark pair if it is energetically favorable. This string breaking continues until all 
particles are connected by strings with sufficiently low tension and are then associated into 
hadrons. In the process of breaking the strings and creating quark pairs, the number of 
particles that are created at small angles with respect to one another increases from that 
which was created in the perturbative parton shower. Thus, hadronization increases particle 
production at small angles, effectively increasing the strength of the collinear singularity and 
decreasing the value of the average angular structure function. 

The effect of hadronization decreasing the average angular structure function can also be 
quantitatively studied. Note that the angular correlation function is just the (squared) mass 
of a jet from constituents that are separated by angular scale R or less. Dasgupta, Magnea 
and Salam [30] studied the effect of non-perturbative physics on the transverse momentum 
and mass distributions of jets at hadron colliders. For the mass, they found that the leading 
correction due to hadronization is 



where C is independent of Rq, the jet radius. For the angular correlation function, we 
expect that the effect of hadronization would also result in a correction proportional to R, 
the resolution parameter of the angular correlation function. We can write 



where Cpert is the perturbative contribution to the angular correlation function and Cnon-pert 
is the non-perturbative contribution. The average angular structure function that follows 
from this is 




(55) 



non-pert 



R, 



(56) 
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< a , 



(57) 
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where the inequahty follows when a > 1. Note that the perturbative angular structure 
function is approximately a and so, indeed, hadronization effects decrease the value of the 
angular structure function. 

The argument presented here and in |30] relies on the one-gluon approximation to de- 
termine the effect of hadronization. Universality of the hadronization and power correction 
effects was argued with the one-gluon approximation in refs. [311 122] and demonstrated for 
event shapes in SCET in refs. [331 [31] • The arguments in refs. [331 El] relied on the boost 
invariance of the soft function for back-to-back jets. How the argument might extend to an 
arbitrary number of jets in arbitrary directions is unclear as the boost invariance is, at least 
naively, broken. We will not discuss how this might be extended, but we note that, because 
of the qualitative and quantitive arguments from the one-gluon approximation, we expect 
that the universality holds in SCET. 



G. The Angular Correlation Function at the LHC 

Finally, we will discuss how the results obtained here for the SCET resummation might 
be extended to the LHC, to processes initiated by pp collisions. For an observable O that 
factorizes at hadron colliders, the cross section can be written in the schematic form 
da 



S{0;f,) (58) 



The beam functions Bi encode the properties of the initial parton i and the matrix Cab 
weights the colliding partons by the appropriate cross section. Indices a and b are implicitly 
summed over. In this case, the flavor of the jet functions depends on the flavor of the initial 
colliding partons which affects the admixture of quark and gluon jets that contribute to 
O. Note also that the soft function includes contributions from radiation from initial state 



partons. Therefore, while not necessarily manifest in Eq. [58| the beam functions implicitly 
affect the jet and soft functions. 

Nevertheless, we expect that the angular correlation function has nice factorization prop- 
erties at hadron colliders. With the goal of computing the average angular structure function, 
we can again ignore anything in the factorization of the angular correlation function that is 
independent of Qa or the resolution parameter R: 



dg, 
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where we have chosen to measure Qa in jet i in an event with N jets. Dependence on the 
beam functions has been retained, however. This is because, for a given set of jets 1, . . . , A'^, 
different initial states contribute to the cross section with different weights. Thus, the beam 
function contribution to the factorization, CabBa{^j)Bh{^)^ is actually not an overall constant 
factor and so must be included. Note also that the color of the colliding partons affects the 
radiation included in the soft function. The beam functions are universal and so can be 
computed once and for all. While this is not a rigorous proof of factorization of the angular 
correlation function, many of the results obtained in the e^e" collider context should be 
able to be recycled for the hadron collider case. This deserves significant future study. 



III. COMPARISON TO FIXED-ORDER CALCULATION 

In this and the following section, we will focus most of our attention on the (proper) 
angular correlation function with a — 2: 

g,(R) = g{R) = ^ J] E,E^ sin % tan ^Q{R - e,j) . (60) 

To evaluate this jet observable in SCET, we must choose the hard, jet and soft scales. For 
many of the comparison plots we choose the following scales: 

Hh^^ , fJ'j^ ^G'^''^ , fJ'S = • (61) 

tan ^ 

These choices of scales minimize logarithms that appear in the resummed distribution. How- 
ever, it is important to understand the dependence of the result on the choice of these scales 
and so we will also present plots in which the scales are varied by the standard factors of 
2 and 1/2. The evaluation of the average angular structure function from the SCET cross 
section is done numerically. Note that for consistency of the factorization the jet scale /ij 
must be larger than the soft scale /is which is the requirement that 

tan|»e;V2. (62) 

To maintain this separation, the resolution parameter cannot be too small; we will only 
consider R > 0.1. For smaller values of R, logarithms of it! become large and must be 
resummed, which is beyond the scope of this paper. 
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FIG. 4. Plots of the average angular structure function for quark (red) and gluon (blue) jets. 
Fig. |4(a) compares the curves from SCET resummation (solid) to anti-fc-r jets from PythiaS 
(dashed). The PythiaS curves were computed from 3 jet final states in which all jets had equal 
energy. Figs. [4(b)|[4(^ and [4(d) compare the PythiaS curves to SCET bands in which the hard, 
jet and soft scales have been varied by a factor of 2. To make these curves, the jet radius has been 
set to be i?o = 1.0 and the energy of the jets is 300 GeV. 



The average angular structure function as computed in SCET is plotted in Fig. |4] where 
the curves for quark and gluon jets are compared to the output of PythiaS. The PythiaS 



curves will be discussed in the next section. Fig. 4(a) compares the quark and gluon curves 



with the hard, jet and soft scales set to their values in Eq. [61] The observations from 
the previous section are apparent with the quark average angular structure function less 
than the gluon average angular structure function and both slightly less than 2. The scale 
variations of these curves are shown in Figs. |4(b) 4(c) and 4(d) Note that in particular 
there is relatively wide range over which the angular structure function varies when the jet 
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and soft scales are changed by a factor of 2. 

Because the average angular correlation function is defined by integrating over the entire 
range of Qa, its value and shape is sensitive to radiation in all regions of phase space. 
Resummation is necessary for an accurate description of the physics in the singular regions 
of phase space while higher fixed-order contributions are necessary for a good description 
in the non-singular regions of phase space. A proper treatment of resummation and fixed- 
order involves consistently matching the two contributions so that the resulting distribution 
is accurate order by order in as over the entire phase space. This matching is a non-trivial 
procedure and, instead, we will just focus on the contribution from higher fixed-order matrix 
elements. This will give us a sense, at least, for how fixed order and resummation affect the 
average angular structure function. 

To do this, we use NLOJet++ v. 4.1.3 [HHl EH], based on the dipole subtraction method 
of 110], to compute the average angular structure function to NLO in e~^e~ collisions. NLO- 
Jet++ can compute matrix elements to NLO for up to 4 final state partons (and, at tree 
level, up to 5 final state partons) and so, by demanding jet requirements, produces jets with 
very few constituents in them. This results in very inefficient calculation of cross sections. 
Also, the public version of NLOJet++ does not record flavor information of partons so the 
identity of quark and gluon jets cannot be easily determined. Further, it is not enough 
that the cross section differential in the angular correlation function at fixed R is smooth 
for the average angular structure function to be smooth. The distributions must also be 
smooth over R so that the derivative that defines the average angular structure function is 
well-behaved. To assuage these issues, in this section, we will define an event-wide angular 
correlation function, where the sum in Eq. [5] runs over all particles in the event. 

The event-wide angular correlation function is defined over all particles in the event with 
no jet algorithm cut. In the limit that there are three final state particles this reduces 
precisely to the angular correlation function of the hardest jet, extending up to an R of 
about the radius of the hardest jet. The angular correlation function will only include the 
contribution from the two closest partons because the third parton must be very far away 
in angle. This argument doesn't hold at higher orders, but for those cases we expect that 
the event-wide definition will be an average over the angular correlation functions of quark 
and gluon jets. 

We present the calculation of the average angular structure function from NLOJet-|--f 
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FIG. 5. Comparison of NLOJet++ calculation of the event-wide average angular structure function 
in 3 jet final states to NLO (dotted) to SCET NLL resummation of the average angular structure 
function for quark jets (solid). Two curves from PythiaS are shown: the dashed curve is the 
average angular structure function for quark jets from e~^e~ — )• 3 jets and the dot-dashed curve is 
the average angular structure function from e^e~ — )• 2 jets. 

for three final state partons to NLO in Fig. [5j The center of mass energy is taken to be 600 
GeV in e^e^ collisions. At a e^e^ collider, most of the time, the hardest jet will contain a 
quark and a radiated gluon so we compare the output of NLOJet++ to NLL resummation 
results for quark jets. Fig. |5]also contains two curves of quark jets from PythiaS which will 
be discussed in the next section. The NLO calculation of the angular structure function is 
approximately fiat and greater than 2 which we interpret as an effective weakening of the 
collinear singularity due to the presence of wide-angle radiation. The fact that the NLO 
result is slightly larger that 2 was anticipated in where it was shown that a jet with 
three constituents should have an average angular structure function larger than 2 by a term 
proportinal to a^. Matching the calculations from NLL and NLO would produce a curve 
that interpolates between the NLL result at small R and the NLO result at large R. 

To generate the NLOJet++ curve, about one trillion events were processed over about 
1 CPU year. Even with this many events, the average angular structure function from 
NLOJet++ is still quite noisy. However, the noise can be reduced by averaging the angular 
structure function over a small range in R at each point. This was done for the curve in 
Fig. |5| Computing the average angular structure function in 4 jet final states to NLO was 
attempted in the same CPU time as the 3 jet results. However, the resulting curves were 
much too noisy to be used. To produce curves at higher orders using NLOJet++ probably 
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requires centuries of CPU time for distributions to converge. However, other programs such 
as BlackHat [H] might be better-suited to higher multiphcity final states at NLO. Work in 
this direction is ongoing. 

IV. PARTON SHOWER MONTE CARLO COMPARISON 

In this section, we compare our calculation of the average angular structure function 
from SCET to the output of Monte Carlo event generator and parton shower. Through the 
Sudakov factor which dictates the probability that no branchings occur between two scales 
of an evolution variable, the parton shower resums logarithms of the evolution variable 
that arise from soft and collinear emissions. Monte Carlo generators create fully exclusive 
events and so the process of resummation of the logarithms is distinct from that in SCET, 
for example, and examining the differences is interesting. The parameter that defines the 
evolution in the parton shower is also (relatively) arbitrary and different choices of the 
evolution variable lead to different emphases on soft or collinear splittings. In addition, 
hadronization and other non-perturbative physics is described by phenomenological models 
which can be used to understand the size and effect of power-suppressed contributions to 
observables. All of these points and their effects will explored in this section. 

For most of the Monte Carlo comparison, we first generated tree-level events for the 
process e^e~ — )■ qqg using MadGraphS v. 1.4.5 |12] at center-of-mass energy of 900 GeV. 
These partons are required to each have equal energy E = 300 GeV so that they are well- 
separated and factorization-breaking terms in the soft function are minimized. These events 
were then showered using the pT-ordered shower of PythiaS v. 8.162 [13]. All default settings 
of Pythia8 were used except for turning hadronization on and off to study the difference. 
In most plots hadronization in Pythia8 has been turned off. To study the effect of using 
different evolution variables in the parton shower we shower the MadGraph events with 
VINCIA V. 1.0.28 IH]. From the showered events, jets were found with the Fast Jet v. 3.0.2 
[15] implementation of the anti-^T algorithm |46j. We choose the jet radius to be Rq = 1.0. 
The three hardest jets are required to have energy between 250 and 350 GeV and we identify 
jets as coming from a quark or gluon by demanding that the cosine of the angle between 
the jet axis and the direction of a parton from MadGraph be greater than 0.9. 

In Fig.|4} we plot the average angular structure function for quark and gluon jets identified 
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in PythiaS (with no hadronization) and the angular structure function as computed in SCET. 
Note that the average angular structure function as computed in SCET is significantly 
smaller than that from PythiaS, especially at larger R. This difference can be attributed to 
higher order effects which were shown in the previous section to increase the value of the 
average angular structure function. Fig. |4] also illustrates the distinction between quark and 
gluon jets. For most of the range of < i? < 1, the average angular structure function for 
gluon jets is greater than that for quark jets, reflecting the fact that gluons have more color 
and radiate more at wider angles than do quarks. This effect is present in both the PythiaS 
curves and the resummed calculation. Because the SCET calculation only included effects 
from jets with at most two constituents, the curves terminate precisely at the jet radius of 
Rq = 1.0. For these anti-fcy jets in PythiaS, the edge effects from the jet algorithm are small, 
extending only over a range of at most R = O.S to R = 1.2. Also, we have not plotted the 
SCET curves below R = 0.1, where they begin to deviate substantially from their value at 
larger R. 

Fig. [5] compares the average angular structure function from NLOJet++ to quark jets 
in SCET and two different curves from PythiaS. The different PythiaS curves exhibit the 
affect of wide angle radiation captured by the jet on the average angular structure function. 
In that figure, the dashed curve is the quark jet average angular structure function from the 
PythiaS sample described above. The dot-dashed curve is the the average angular structure 
function from quark jets from e^e^ — > qq samples generated and showered in (otherwise 
default) PythiaS. The center of mass energy was set to be 600 GeV and the jets were 
required to have energy within 50 GeV of 300 GeV. Higher order effects are obvious. The 
PythiaS curves agree well with one another at small i? up to i? ~ 0.6 and then diverge at 
larger R. The jets in the 3 jet sample collect wide-angle radiation from the neighboring jets 
which increases the average angular structure function at large R. To fully understand the 
rise within an analytic calculation requires matching fixed-order to resummed result. Fixed- 
order contributions are responsible for the wide-angle emissions that increase the average 
angular structure function because SCET factorization effectively decouples the jets. 



As discussed in Sec. II F we expect the effect from non-perturbative physics on the 
angular structure function to be small and relatively well-understood. In particular, relying 
on arguments from the one-gluon approximation, we expect that hadronization increases the 
strength of the coUinear singularity and that this effect is most prominent at small values of 
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FIG. 6. Comparison of the average angular structure function as computed in PythiaS with (dotted) 
and without (dashed) hadronization. 

R. In Fig. |6| we have plotted the average angular structure function for quarks and gluons 
comparing the curves with hadronization turned on or off. Indeed, the effect is small but 
unambiguous: hadronization effectively increases the strength of the coUinear singularity. 
As discussed earlier, extending the arguments from [33| IM] on the effect of non-perturbative 
physics would be greatly desired to fully describe (at least) average behavior of hadronization 
for events with multiple jets. 



A. Monte Carlo Error Estimates 

Finally in this section, we would like to get a handle on the error or uncertainty in 
the Monte Carlo parton shower in Pythia. Typically, this is done by studying the output 
of different tunes of the same Monte Carlo program or comparing different Monte Carlo 
programs altogether. In particular, as is relevant for the parton shower, the evolution 
variable of the parton shower dictates when and how emissions should occur. Ref. [17] 
observed differences in event shape variables as computed in Pythia 6.4 [JS] between two 
tunes; one py-ordered and the other virtuality ordered. However, these two tunes had 
other distinctions as well and so purely the effect of the evolution variable is obscured. Also, 
comparing two different Monte Carlos is subtle because the number of differences is typically 
huge and so isolating effects of single parameters or choices is very difficult. 

Here, we would like to study the effect of different evolution variables in the parton 
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FIG. 7. Comparison of SCET computation (solid) and PythiaS (dashed) of average angular struc- 
ture function to the output of VINCIA Monte Carlo parton shower with two different evolution 
variables: pT and virtuality. The shaded region lies between the curves from VINCIA. 



shower. The choice of the evolution variable is only a change of variables in the Sudakov 
form factor and so must produce the exact same leading-log resummation for any (consistent) 
choice of evolution variable. However, the choice of evolution variable can lead to higher 
log-order effects through the scale at which is evaluated or by emphasizing soft over 
coUinear splittings, for example. To study the differences, we use the VINCIA [44J parton 
shower plug-in for PythiaS which is based on 2-to-3 splittings as opposed to the standard 
l-to-2 splittings as in Pythia and Herwig |19] . VINCIA includes a flag which allows the user 
to change only the evolution variable. For concreteness, we will consider p^-ordering and 
virtuality ordering. 

In Fig. [7] we have plotted the SCET resummation and PythiaS output for the average 
angular structure function as well as a band which extends over the range between the 
output of the pT-ordered and the virtuality ordered shower in VINCIA. The exact same 
requirements on the jets were made in the VINCIA sample as in the PythiaS sample as 
described earlier. Over most of the range in R, the lower edge of the VINCIA band is set 
by the p^-ordered shower and upper edge by the virtuality ordered shower. This is expected 
as the pT-ordered shower emphasizes collinear emissions more than the virtuality-ordered 
shower. Note also that the band is slightly above the output of the pr ordered shower in 
PythiaS. Part of this effect could be due to the default matrix element matching in VINCIA: 
final states with up to 5 partons are matched to tree-level matrix elements. Regardless of 
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the details, the effect of changing the evolution variable is large. Understanding if and how 
parton showers resum higher order logarithms with different evolution variables, matching 
schemes, etc., is necessary to understand the source of the differences. 



V. CONCLUSIONS 

The average angular correlation and structure functions capture the average scaling prop- 
erties of QCD jets. We have presented a calculation of the angular correlation function to 
NLL accuracy in SCET and compared this result to the PythiaS Monte Carlo parton shower 
and to fixed-order results from NLOJet++. Comparing the resummed SCET result to the 
fixed-order NLO Jet++ result provides a good understanding as to the behavior of the parton 
shower result. However, for a full understanding, matching of the resummed and fixed-order 
distributions is required. Much like the jet shape [15], the average angular structure function 
could be used for tuning of the Monte Carlo. Because it is a two-point correlation function, 
the angular correlation function captures distinct information from the jet shape and so this 
tuning would be non-trivial. 

There are several directions for extending the study presented here. First, it would be 
desirable to compute the angular correlation function in collisions at the LHC. It remains an 
outstanding problem to use SCET to resum logarithms for arbitrary observables in hadron 
colliders because factorization of the (colored) initial and final states is highly non-trivial. 



However, using the observations from Sec. II G the computation of the average angular 



structure function at hadron colliders might only require a reinterpretation of the results 
presented here. Recently, NLO results were obtained for pp — )■ 4j events [50j from which 
any IRC-safe observable could be computed. In particular, for four final state partons, 
the hardest jet can contain up to three constituents which would be beyond the resummed 
order in the SCET calculation. Also, at a hadron collider, underlying event or pile-up 
produce significant background radiation that can be collected into a jet. A procedure 
to determine the contribution to a jet from these non-perturbative sources is necessary to 
properly determine jet energy scales and to study substructure. The results presented here 
could be used to determine the average contribution to a jet using the procedure introduced 
in [ig. 

For a more accurate prediction of the angular correlation function, matching of NLL and 
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NLO results must be done to have good control of the distribution over the entire phase 
space. Factorization of jet observables allows for a process-independent computation of 
the NLL resummed result; however, the fixed-order calculation is process dependent and 
must be couched in a particular study. We showed that the average angular structure 
function is sensitive to wide-angle radiation so matching is vital for accurate predictions. 
NLOJet++ or results like those from [50] are promising in their applicability to generic 
processes. It is unlikely that QCD jet observables can be reliably computed to NNLL or 
beyond analytically because non-global logarithms become important. Nevertheless, by 
studying limiting behavior such as in |B] the effect of these non-global logarithms might be 
reduced. 

Finally, as there exist few jet substructure observables that have been (or even can be) 
computed with analytic methods, it is important to compute those that are possible. The 
calculation of the angular correlation function provides powerful insight into the behavior 
of QCD and the dynamic properties of jets. Though scale-invariance is broken in QCD by a 
running coupling, jets maintain a fractal, conformal structure to very good approximation 
over a wide dynamical range. 



Appendix A: Measured Jet Functions 

Here, we present the finite pieces of the measured jet functions for quark and gluon jets 
as defined by a kx-type algorithm for the angular correlation function. These functions are 
composed of contributions from 5-functions and +-distributions. For a function g{x), we 
define the -I— distribution as [51] 

[g{x)e{x)]^ = g{x)e{x) - S{x) f dx' g{x') , (Al) 

so that 

/ dx [g{x)Q{x)]^ = . (A2) 
Jo 

From this definition, it is straightforward to compute the measured jet functions. The terms 



that are infinite in four dimensions were presented in Sec. [11 B[ The terms that are finite in 
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four dimensions are, for a quark jet: 
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is the largest value that Ga can take for a jet with two constituents: 
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where we have taken the leading A dependence for the coUinear modes. 



Appendix B: Resummed Distribution for Angular Correlation Function 



The expression for the resummed cross section is 
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wj, 00 s, Kj and Ks result from the resummation of the individual jet and soft functions 
[n2HSS]- The functions uj and ujk are defined by ojj = up{fi,fij) and us = —ujf{^,^s) 
where 
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/3o/ 47r 
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(B2) 



/3o is the coefficient of the one-loop /3-function as defined in Eq. 30 and (3i is the two-loop 
coeffcient: 



(B3) 



r is the ratio of the strong coupling at two scales: 



tts(/io) 

and the energy dependence of the strong coupling is given by the two-loop expression 
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for as evaluated at the two scales and Q. The terms r[!^gp and Fj^^p are the one- and 
two-loop coefficients of the cusp anomalous dimension. Their ratio is given by [57] 
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where Cj is defined in Eq. 49 Ks is defined similarly 
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The functions fj and fs are generated by the convolution of the jet and soft functions. 
Accurate to NLL, they are 
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is the harmonic number function defined by 
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and tjj^^^x) is the trigamma function 
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Note that the logarithms in these functions can be minimized by choosing 
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For expansion of the resummed distribution, the following relations are needed: 
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